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Abstract 

We compute the 0(/5oa^) QCD corrections to the Standard Model decay 
t W~^b as weh as the non standard decay t H^b. We then use our 
resuhs to compute the BLM scale for these decays, and study the dependence 
of the BLM scale on the mass of the boson in the decay. We find that the BLM 
procedure gives extremely small scales when top decays into heavy bosons. 
When we try to improve the BLM scale by reexpressing rates as a functions 
of the top quark running mass, we find the BLM scale exhibits unnatural 
dependence on the boson mass, casting doubts on the applicability of the 
BLM procedure to these processes. 
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The top quark, recently observed by CDF and DO 0, will be the focus of much 
experimental attention in coming years. It will be important to test the Standard Model 
predictions for the width of the top quark as well as search for rare or non-standard decays 
which may give us a glimpse into post-Standard Model physics. Because of this it is im- 
portant to have precise theoretical predictions for top quark decay, both for the Standard 
Model and for theories which propose to go beyond it. 

In this paper, we calculate the Oi^Poa"^) corrections to the decays t W~^b and t —>■ H~^b, 
including the full dependence on the mass of the boson in the decay. (We will neglect the 
mass of the b quark, however.) In many QCD calculations, these corrections form the 
dominant part of the complete 0(a;^) correction. As examples we present the perturbative 
corrections to R{e^e^ — » hadrons), r(r Vr hadrons) and the series relating the pole mass 
and the MS mass of a heavy quark: 
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In each case the 0(/?o«s) term is a good approximation to the complete 0{a1) correction. 
Because of this, the calculation of the 0{j3oa'^) corrections for top quark decays is of interest 
in its own right. This calculation is a principal result of this paper. 

The 0{Poa'^) corrections can also be used to set the scale of the QCD coupling, via 
the scheme of Brodsky, Lepage, and Mackenzie Varying the scale of the QCD coupling 
changes the n*'^ order coefficient in the perturbation series by a term proportional to /3q 
Since (3q is numerically rather large, it is clear that a poor choice in scale can result in a 
poorly behaved perturbative expansion. The BLM scheme sets the scale so that the 0(/5oa^) 
term is cancelled. 

The BLM scale is commonly used as a criterion for assessing the validity of a perturbative 
QCD series. A small BLM scale is often interpreted as a sign that higher order corrections 
are not well behaved and that the perturbation series is not to be trusted. However, for the 
top quark decays discussed in this paper, we will see that the BLM scale can become much 
smaller than any of the relevant mass scales in the problem. For instance, we will find a 
BLM scale of 200 MeV for a 175 GeV top quark decaying into 150 GeV Higgs and massless 
b quark. Since the energy released in the decay is far greater than the QCD scale, it is hard 
to believe that perturbative QCD is really in trouble. Indeed, by examining the coefficients 
in the perturbative series we will see that the small BLM scale is a consequence of a small 
leading order coefficient rather than a large second order coefficient. For top quark decaying 
into very heavy Higgs, the BLM scale ceases to be a useful tool in analyzing the validity of 
the perturbative QCD series. 

More sophisticated attempts at getting at the "true" scale of a leading order calculation 
were made in refs. These authors resum renormalon chains, which include the entire 
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/9o~^q;" series, for processes such as R(e"'"e~ hadrons), r(r Ur hadrons), and r(b — >• 
cev). This can be viewed as an all orders generalization of the BLM scheme, which only 
sums the 0{l3oa1) term correctly. The sum of this infinite series, defined as the principal 
value of the Borel integral, can be absorbed into the scale of the coupling constant in the 
leading order calculation. The efffective scale found in these calculations rarely differs from 
the BLM scale by more than ±50%, indicating that the BLM scale may be used to obtain 
a rough estimate of the sum of this subclass of higher order corrections. The quality of the 
perturbative series is not measured by this effective scale, however, but by the imaginary 
part of the Borel integral, which is interpreted as being the intrinsic uncertainty which one 
must encounter when trying to sum a series which is not Borel summable. In ref. 0, a 
renormalon summation is done for the decay t W^h in the limit My/ —>■ 0. The authors 
of find an effective scale of 0.89 rrit, to be compared to the BLM scale 0.122 rrit [0. 
The uncertainty in the summation is estimated to be less than one percent, in line with 
our intuition that strong interaction corrections to this decay should be under perturbative 
control. It will be of interest to redo the calculation allowing for finite W and Higgs masses. 



Such a calculation is currently in progress [12 



We begin by considering the Standard Model top quark decay t bW^. The QCD 
corrections to the tree level decay rate have been computed to 0{as) in ref. 0: 



2a, (/i) / 2tt^ 5 



T{t W+b) = To 
where xw = ^w/^h ^o, the tree level decay rate, is 



(4) 



To = I'f^ (1 + "^^w) (1 - xwf = 1.33 GeV. (5) 

In eqns. (^) and (^, we have neglected the mass of the b-quark. We use nit = 175 GeV, 
mw = 80 GeV, and Vtt = 1. For these values of mw and rrit, xw = 0.2; 5(0.2) = 0.31. 

Smith and Voloshin have derived the following integral representation for the 0{(3oO!.'^) 
term in top quark decay: 

^ Joa^'Km) f- _ _^r^r,^o)] (6) 

Here, r'^^)(yu) is the one loop correction computed as if the gluon had mass /i and ^^^•'(mj) is 
the strong coupling constant defined in the V-scheme of Brodsky, Lepage, and MacKenzie . 
The V-scheme coupling constant is related to the MS coupling constant by 

«i''H/i) = «Pne-M(l + OK)). (7) 
Computing in the limit of my/ —>■ 0, Smith and Voloshin find 

5r(^) = -1.73^0 (^^)'. (8) 

Their result is confirmed in ref. |^ , where the calculation is done by direct evaluation of the 
Uf dependent 0(a;^) diagrams. Smith and Voloshin use this result to find a BLM scale of 
0.122 rrit (in the MS scheme) for top quark decay. 
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FIG. 1. Coefficients in perturbative expansion for t W^h; f{xw) - solid, g{xw) - dashed. 



In this paper, we redo the calculation of 0(/?o«s) corrections to t — > W'^h including the 
effect of finite W boson mass. We calculate the 0{as) corrections to t — W^h, giving the 
gluon a finite mass. Integrals over Feynman parameters in the virtual gluon correction and 
over phase space in the brehmstrahlung graphs are computed numerically. The infrared 
divergences in each of these calculations cancel in their sum. Taking the gluon mass to zero 
we find our numerical result agrees with the analytic expressions in ref. [0. We then use 
eq. (H) to compute the 0(/3oag) corrections to t — > W^h. The rate is given by 



T{t W+h) = To 



1 + f{xw) ^ 9 [xw) Po- 



tt 



(9) 



where the functions f{xw) and g{xw) are shown in Fig. |l]. The analytic expression for 
f{xw) can be found in ref. 0. For xw = 0.2, f{xw) = —2.51, g{xw) = —1.96. 
The BLM scale, given by 



^■blm{xh) = m exp 



'^9{xh) 



(10) 



is plotted in Fig. |^ as a function of xw For xw = 0.2, we find fisLM = 0.09 m^. If we use 
this scale for in eq. @), we find that QCD corrections reduce the width of the top quark 
by approximately 13%. (For Og we use the one-loop expression for the running coupling 
constant in the MS scheme, setting nf = 5 and Aj-fg = 110 MeV.) If we choose to leave the 
scale for the coupling equal to rrit and include the O^Poa^) term, the top quark decay rate 
is reduced by 11%. The effect of BLM scale setting on the theoretical prediction is small 
because for small xw, ctgifJ^BLAi)/'^ is still perturbative. 

Our computation of the BLM scale agrees with a similar calculation in ref. 0, which 
studies the perturbative QCD corrections to the semileptonic decay of the bottom quark. 
These authors compute the BLM scale for the differential rate dT/dq"^, where is the 
invariant mass squared of the lepton pair. Since the QCD corrections to the semileptonic 
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FIG. 2. BLM scale for top quark decays. The dotted line is for the decay t W^b; 
X = m^^/mf. For m^^/mf = 0.2, fiBLM = 0.09 nit. The solid line is for t H~^b; x = mj^/mf. 



decay are identical to those for the decay of top into a b quark and real W boson, we can 
simply substitute rrit for rub and m^r for q"^ and recover their result from ours. 

Next we would like to consider the QCD corrections to the decay of the top quark into 
a b quark and charged Higgs boson. A charged Higgs boson appears in any extension of 
the Standard Model which has more than one Higgs doublet. In general, it can have scalar 
or pseudoscalar couplings, so we parametrize the Feynman rule for the t — — b vertex 
as iu{h){s + p75)M(t), where the parameters s, p depend on the specific model for the Higgs 
sector. The tree level decay rate (in the limit 0) is 



Vo{t ^ H+h) 
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The 0{as) corrections to this decay were calculated in ref. [|T^]. 



Tit ^ H^) = Fo 



1-Xh Xh 



(12) 



-5 ln(l -xh)-2 U2{xh) + 2 Li2(l -xh) + -- vr^jj , 

where xh = m'jj/mf. The procedure for computing the 0(/3o«s) corrections to t — H^b is 
identical to that for t W^h discussed earlier. Our numerical calculations for the one loop 
corrections agree with eq. ([12D when we take the gluon mass to zero. For t H^h decay, 
we define coeffiecient functions /{xh) and g^xu) exactly as we defined f{xw) and g{xw) 
for t — > W^h. /{xh) and g^xu) are plotted in Fig. ||. (The logarithmic divergence of these 
functions as xh ^ 1 is an artifact of the approximation mb = 0.) 

The BLM scale for t H^h is plotted in Fig. ^ as a function of xh- Because the energy 
available for the gluon decreases as the mass of the boson increases, we expect the relevant 
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scale for the QCD coupling to decrease as mw, — mt. What is surprising is that for both 
decays the BLM scale effectively vanishes well before the kinematic limit. For example, if we 
were to use the criterion as^HBLu) /t^ < 1 for a well behaved perturbation series, we would 
conclude that the perturbative computation of QCD corrections to top decay into charged 
Higgs becomes unreliable for xh > 0.75 or itlh > 150 GeV, 25GeV below the kinematic 
limit irit — ~ 175 GeV. 

Upon examining Fig. ^, it is clear that the small BLM scale obtained as — 1 is as 
much a consequence of the first order correction becoming small as it is a consequence of 
second order correction becoming large. It is also useful to look at the the numerical values 
of each term in the series. Considering the case xh = 0.75 we find 



r = rn 



1.22 



31.8 



TT 



OlsWIt) 



To [1 - 0.043 - 0.039] . 



(13) 



The second order correction is almost as large as the first; this is why we get the very low 
BLM scale 200 MeV. Using this BLM scale in the leading order expression we find the QCD 
corrections to the rate are —56%. On the other hand, the absolute size of each correction is 
small, and the expansion parameter a{mt)/Ti is only ~ 0.003. To get a mere 4% correction 
to the total rate would require a third order term of the size ~ 900 a^/vr^ or ~ 15 /^Qa^/vr^. 
This is not unreasonable, but it does require a large coefficient in the third order term. To 
get the kind of correction implied by the BLM method implies huge coefficients at higher 
orders. It seems likely that in this case the BLM scale setting method is overestimating the 
effect of higher order corrections. 

Next, we will attempt to improve the BLM scale by rewriting the rate as a function of 
the top quark running mass. It is shown in ref. [O that the leading renormalon singularity 



in the perturbative series for heavy quark decay is cancelled when the rate is expressed 
in terms of the running mass instead of the pole mass. Since the asymptotic growth of 
coefficients in the perturbative series is controlled by renormalon singularites, one expects 



6 




FIG. 4. BLM scale for t H'^h when rate is expressed as a function of the top quark MS 
mass instead of the pole mass. 



smaller coefficients in the perturbative series if we rewrite the series for top quark decay 
in terms of running top quark mass instead of the pole mass. While it is not obvious that 
arguments based on the asymptotic behavior of the perturbative series are relevant to the 
lowest order terms, it is in fact true that the coefficients of low order 0(/3q a""'"^) terms in the 
decay t W^h (in the limit My/ 0) are significantly reduced when the series is expressed 
as a function of the running mass . 

In the limit — *• 0, the parameters s, p are proportional to m^, so the total rate is 
proportional to m^. Making the substitution in eq. (|l]), then applying the BLM procedure 
to the new perturbative series results in the BLM scale shown in Fig. ^. The BLM scales 
computed for the two series are drastically different. The series obtained when the rate is 
expressed as a function of the MS mass has a higher BLM scale for almost all values of xh- 
However, the dependence of the new BLM scale on xh seems unphysical. The BLM scale 
shown in Fig. ^ actually becomes larger as xh 1- 

To understand why the BLM scale changes so dramatically, we examine the coefficients 
of the one and two loop corrections for the new series. When expressed as a function of the 
running top quark mass, the top quark width is given by: 



71^ 



(14) 



with a BLM scale 



fJ-BLM = mt exp 



2g{xH) +9.36' 
/(xh)+4 



The functions /{xh) + 4 and g^xu) + 4.68 are also plotted in Fig. |^. In this case the ratio 
{gi^xu) + 4.68)/(/(x//) + 4) is small, especially for large xh, and we find much larger BLM 
scales. As xh — > 1, the replacement of the pole mass with the running mass leads to a 
nearly exact cancellation of the two loop correction. Since there is no physical reason to 
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expect smaller higher order corrections in the limit xh 1 it is doubtful that we should 
regard this cancellation as anything other than fortuitous. We see again that the physical 
significance of the BLM scale is questionable. 

In summary, we have computed the 0(/3oQ;^) corrections to top quark decay into W and 
Higgs bosons, and apphed the BLM scale setting method. For the case of decay into heavy 
bosons, we found extremely small BLM scales for processes that we intuitively expect to be 
calculable in pcrturbative QCD. The small BLM scales can be made larger by reexpressing 
decay rates as functions of the top quark running mass rather than the pole mass. However, 
in this case the BLM scale exhibits unphysical dependence on the mass of the boson. These 
problems lead us to believe that for the BLM scale is not a useful tool for analyzing QCD 
corrections to these decays. A more detailed probe of higher order corrections is needed in 
order to set the scale. 

It is a pleasure to acknowledge helpful discussions with Adam Falk, George Chiladze, 
and Konstantin Matchev. This work was supported by the National Science Foundation 
under Grant No. PHY-9404057. 
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